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ABSTRACT 

We use the grid of hydrodynamic accretion disc calculations of Stchle (1999) to 
construct orbital phase dependent emission line profiles of thin discs carrying spiral 
density waves. The observational signatures of spiral waves are explored to establish 
the feasibility of detecting spiral waves in cataclysmic variable discs using prominent 
emission lines in the visible range of the spectrum. For high Mach number accretion 
discs (M = v^/cs ~ 15 - 30), we find that the spiral shock arms are so tightly wound 
that they leave few obvious fingerprints in the emission lines. Only a minor variation 
of the double peak separation in the line profile at a level of ~8% is produced. For 
accretion discs in outburst (Af ci 5 - 20) however, the lines are dominated by the 
emission from an m—2 spiral pattern in the disc. We show that reliable Doppler 
tomograms of spiral shock patterns can be reconstructed provided that a signal to 
noise of at least 15, a wavelength resolution of ~ 80 km s~^ and a time resolution of 
'^ 50 spectra per binary orbit is achieved. We confirm that the observed spiral pattern 
in the disc of IP Pegasi can be reproduced by tidal density waves in the accretion disc 
and demands the presence of a large, hot disc, at least in the early outburst stages. 
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1 INTRODUCTION 

Cataclysmic variables (CVs) harbour a geometrically thin 
viscous accretion disc through which mass from a late type 
companion star is accreted onto a white dwarf. Rapid bright- 
ness variations of several magnitudes displayed by the sub 
class of dwarf novae indicate that information is transported 
along the accretion disc on a viscous time scale of a few days 
(see Warner 1995 for an extensive review). The dynamical 
time scale of disc material on the other hand, is given by 
the local Kepler velocity and ranges from seconds near the 
white dwarf, to a considerable fraction of the orbital period 
(hours) at the outer disc. The effective viscous transport 
displayed by these discs demands the presence of a highly 
efficient source of angular momentum transport within the 
discs, whose physical origin is poorly understood. 

The convenient time scales in CVs allow us to study the 
changing properties of viscous accretion discs in real-time. 
As the accretion disc is viewed from different angles during 
a binary orbit, image reconstruction techniques (IRTs) such 
as eclipse mapping (Home 1985) or Doppler tomography 
(Marsh & Home 1988) are able to provide spatially resolved 
images of accretion discs on a micro-arcsecond scale. Since 



IRTs use the luminosity variation of the disc with binary 
phase during a considerable fraction of the binary orbit, the 
reconstructed images are time averaged pictures of the disc 
pattern. This implies that IRTs can reconstruct those disc 
features which are stable over a period comparable to a dy- 
namical time scale in the co-rotating frame of the binary. 

Numerical calculations by various authors (e.g. Sawada, 
Matsuda & Hachisu 1986, Rozyczka & Spruit 1993, 
Heemskerk 1994, Larwood 1997, Yukawa, Boffin & Matsuda 
1997, Stehle 1999) show that tidally induced waves can re- 
sult in two prominent spiral shock arms, which are station- 
ary on a dynamical time scale in a co-rotating frame. The 
shocks are excited in the outer disc by the tidal action of 
the secondary star. In the binary potential field, the disc 
gas revolves around the white dwarf on non-axisymmetric 
orbits where the streamline-density for a pressure free gas 
is highest towards the secondary and lowest normal to it 
(Paczynski 1977). In a low-temperature disc, the shape of 
the streamlines are still similar to the pressure free gas, but 
the variation of the streamline-density with azimuth along 
a particle orbit will excite pressure waves which steepen to 
shock waves while they travel inward (Spruit, Matsuda, In- 
oue & Sawada 1987). Additionally, shocks can be excited by 
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the crossing of particle orbits if the disc is larger than the 
tidal radius. The opening angle of the shock arms depend 
mainly on the azimuthal Mach number M = v^/cs with u^ 
the azimuthal disc gas velocity and Cs the sound speed, i.e. 
on the rate with which the shock is sheared by differential 
rotation of the disc gas while the shock travels inward radi- 
ally (Rozyczka & Spruit 1989). 

A strong radial viscous force, if present in accretion 
discs, may on the other hand damp the shocks over a short 
range or may even prevent the pressure disturbance to steep 
into a shock firsthand (Savonije, Papaloizou & Lin 1994). 
As the origin of the viscosity is not known, we neglect in 
our calculations any radial viscous forces and thus derive 
a maximum amplitude for the shocks. If this is a sensible 
assumption is one of the main questions to be answered by 
comparing accretion disc models with yet to come detailed 
observations. 

Despite the firm prediction of shocks in accretion discs, 
two armed spiral signatures have only recently been found by 
Steeghs, Harlaftis & Home (1997) and Harlaftis et al. (1999) 
in the eclipsing CV IP Pegasi during outburst. Doppler 
tomograms revealed a strong two-armed spiral pattern in 
the outer disc. Even though Doppler tomography is ide- 
ally suited to detect asymmetric disc structures (see e.g. 
Robinson, Marsh & Smak 1993), it is unclear if previous 
observations failed in reconstructing spiral arms because of 
insufficient data quality or because spiral arms where not 
present in the disc at the time of the observation. To clarify 
this question the purpose of this paper is to predict obser- 
vational signatures of tidal density waves in high M and low 
M accretion discs, the latter appropriate for accretion discs 
during dwarf nova outbursts. 

We use the emission line model of Home (1995) to cal- 
culate emission line profiles for various binary configura- 
tions, using the grid of full hydrodynamic disc calculations 
by Stehle (1999) as underlying models. 

The paper is organized as follows: in Section H we will 
introduce our time averaged model calculations and our 
method of constructing Doppler tomograms and line pro- 
files. In Section H we present model predictions for several 
relevant cases. The effect of instrumental resolution and sig- 
nal to noise on the detectability of tidal structure is explored 



in Section 0. We conclude the paper in Section 
cussion and summarize our results in Section 



with a dis- 



2 EMISSION LINES FROM SPIRAL SHOCKS 

2.1 The grid of hydrodynamic disc models 

As our underlying disc models we use the grid of geomet- 
rically thin 2D-disc calculations by Stehle (1999). The disc 
thickness is explicitly followed by two additional equations 
in a one-zone model (see Stehle & Spruit 1999 for a de- 
tailed description). The spatial and temporal evolution of 
the disc temperature follows tidal and viscous heating, the 
latter in the a-ansatz of Shakura & Sunyaev (1973), as well 
as radiation from the disc surfaces. 

Our time dependent calculations were time averaged 
over at least one binary orbit in the co-rotating frame of the 
binary, i.e. 
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Figure 1. The effective temperature T^ff{r) versus the disc radius 
r/rij-^ for models I (a = 0.01, lower line) and II (o=0.3, top 
line). Ti,{i{r) follows closely the power r~^''^ (dashed line), as 
expected for hot, stationary accretion discs. Near the outer disc 
edge, the effects of tidal heating and non-axisymmetry result in 
a wide range of disc temperatures. 
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where / is any conservative physical quantity, like the mo- 
mentum in radial direction pr, the angular momentum p^, 
the vertical momentum pz, the surface density E or the in- 
ternal energy density e. Other quantities, as e.g. the disc gas 
velocity, are derived from the time averaged values. We will 
deal further on only with time averaged values and thus we 
subsequently drop the bar in our notation for convenience. 
By averaging our models in time, only those disc features 
which are reasonably stable over a dynamical time scale will 
still be visible. 

We will be presenting the results of two model simula- 
tions. Both models are simulations of the accretion disc in a 
binary with mass ratio q — M2/M-1 — 0.3 and orbital period 
P — 2.3 hours. The a-type Shakura & Sunyaev viscosity in 
model I is 0.01 and we time-averaged over ~1.5 binary or- 
bits. Model II represents a typical dwarf nova outburst disc 
with a = 0.3, time averaged over ^ one orbital period. We 
also calculated model predictions for various other mass ra- 
tios and different averaging times which all share the same 
qualitative features since the spiral arms are very stable in 
the binary frame and mainly depend on the disc tempera- 
ture (see also Stehle 1999). These two models in particular 
were selected as a typical comparison between a small, cool 
disc (I) versus a large, hot viscous disc (II) illustrative for 
the discs in CVs. We show in Fig. hi the run of T^g with disc 
radius r in units of the distance to the Li point. In both 
cases Tcff(r-) oc r"^'** in the inner parts of the disc, as ex- 
pected from stationary viscous accretion discs. Close to the 
outer disc edge though, tidal heating and non-axisymmetry 
result in the large spread of Tcff (r). 



2.2 Doppler tomograms and Une profiles 

For our analysis we use a right-handed Cartesian coordi- 
nate system (cx, Sy, e^) centred on the white dwarf, with the 
positive X-axis in the direction of the secondary star and 
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Bz parallel to the rotation vector of the accretion disc. The 
disc gas velocities are transferred to inertial velocities cor- 
responding to the centre of mass of the binary. To calculate 
the line emissivity for each disc grid cell, we follow the lines 
of Home & Marsh (1986) and Home (1995). The emission 
line surface brightness J is obtained by integrating over the 
local line profile times the foreshortening due to the orbital 
inclination i; 



J = coai / 51,(1 — exp (—r^)) d// 



(2) 



In the case of optical thin line emission, J(x, y) is propor- 
tional to the surface density T.{x, y) and the local line source 
function Sl{\)- However, most strong emission lines such 
as the hydrogen Balmer series are optically thick and satu- 
rated. The area under the profile is then well approximated 
by the product of Sl and the frequency interval over which 
the optical depth exceeds unity. Using eq. (15) from Home 
(1995), the local line emissivity of a saturated emission line 
with rest wavelength Ao is in that case given by: 



AV 



J(x,y) = Sl—, — cosi-\/81nro 
Ao 



(3) 



AV is the total velocity dispersion along the line of sight, 
TO the optical depth of the centre of the emission line. If we 
furthermore assume local thermodynamic equilibrium, Sl 
equals the Planck function. The ve locity disp ersion consists 
of a thermal component Kh = yJkT /m-nA, with mn the 
mass of a hydrogen atom and A the atomic weight, and 
a shear term 14hcar, derived from the full, non-Keplerian, 
velocity field. /S.V is evaluated for each grid cell and viewing 
angle according to 



/\v^ = Kh + vL 



(4) 
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e denotes the earth vector pointing from the grid cell towards 
the observer, and VF is the tensor of the local velocity gra- 
dient including all anisotropic terms. In this way the velocity 
field of our simulations self consistently provides the total 
velocity dispersion along each line of sight. AZ denotes the 
vertical extent of the emission line layer and is assumed to 
heO<AZ<H. 
Since tq oc T./AV we derive: 



J{x,y,ct),i) oc B^(roff) AV{x,y,4>,i)J\-a{ 



T.ix.y) 



•■ AVix, y,<j,,i) 



1(7) 



The line emissivity is mainly sensitive to the local velocity 
dispersion AV. The large gradients across the spiral shocks 
will therefor enhance the local emissivity of saturated lines 
significantly at the location of the spiral arms. 

The line emissivity is transformed from position (a::, y) 
to velocity coordinates (14, Vy) by multiplying with the Ja- 
cobian of the coordinate transformation, i.e. 
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For each binary phase (j) the line profile f{v, (f>) is synthe- 
sised by adding up the individual contributions of the grid 
cells where the emission line contribution is Doppler shifted 
according to the local gas velocity along the line of sight, 
Vdop ~ —V-e. We adopt the usual convention that phase 
0.0 corresponds to inferior conjunction of the companion 
star, i.e. mid-eclipse in high inclination systems. The visibil- 
ity of each grid cell is tested for eclipses by the companion 
star which we assume to fill its critical Roche volume. Self- 
shadowing of the disc, important for high-inclination CVs 
due to the vertically extended disc rim will be discussed 
elsewhere (Stehle & Steeghs, MNRAS, in preparation). 

The above description to construct emission line pro- 
files holds principally for all saturated emission lines. We 
choose the Hydrogen /3 line at 486lA as our reference line 
for most calculations. The Balmer lines are the most promi- 
nent lines in CV spectra and allow disc reconstructions at 
a high signal to noise. It serves as a typical representative 
line observed in optical emission line studies of CVs. Since 
we are concerned with the dynamical properties of the line 
profiles rather than calculating absolute line strengths, this 
particular choice does not afi'ect our subsequent analysis. 
To illustrate the effects of saturation we will calculate both 
the case of no shear and maximum shear broadening as well 
as compare with the properties of a heavy element such as 
Calcium. 



3 MODEL PREDICTIONS 

In this section we present predicted emission line profiles and 
Doppler tomograms for the disc models limited only by the 
finite grid resolution of our hydrodynamical calculations. As 
we will see in Sec. q, the limits set by the instrumentation 
will dominate over our limited numerical grid resolution and 
it is thus appropriate to call our tomograms in this Section 
"ideal". We identify the signatures of spiral waves in the 
emission lines of accretion discs. As the global shape of the 
spiral shocks is characterised by the Mach number of the 
disc, we present our model predictions in descending order 
of M. 



3.1 Axisymmetric discs 

The properties of emission lines from axisymmetric accretion 
discs have been investigated by several authors (e.g. Smak 
1981, Home & Marsh 1986), of which we summarise a few 
relevant features. The large Kepler velocities of a few thou- 
sand km s~^ in the inner disc down to several hundreds km 
s~^ in the outer disc, produce very broad emission line pro- 
files, decoupling the global shape of the line from the local 
line radiative transfer. A double peak shape results where 
the peak velocities refiect the motion of the gas at the outer 
disc edge, which in turn can be used to estimate the size of 
the disc if a velocity field, such as Keplerian, is assumed. 

The orbital velocity of the white dwarf is sometimes 
estimated by measuring the centroid of the observed line 
profile, which should show a sinusoidal dependence with or- 
bital phase as the white dwarf orbits around the binary cen- 
ter of mass. While the line profile as a whole thus moves, 
the separation between the double peaks is constant with 
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Figure 2. Model I: an a = 0.01 disc in a binary witli q = M2/M1 = 0.3. Top panels show the distribution of line emission in both spatial 
coordinates (left) as well as in the V^V^-pIane. The crosses in the tomogram denote the projected radial velocity of the white dwarf 
(lower cross) and secondary star. Bottom panels display the calculated emission line profiles for Hydrogen-like lines at an inclination of 
80°. Bottom left with no shear broadening (i.e. a thin emission line layer), right with maximum shear broadening when line emission is 
produced across the whole disc height. To the right of each trailed spectrogram is the total line flux light curve, and below it the emission 
line profile at orbital phase 0.75. 
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Figure 3. The effect of atomic weight on the shear broadened profiles. Left are the line profiles for Model I viewed at an inclination of 
i = 45° for a Hydrogen line (A=l). Right panel is for a saturated Calcium line (A=40). Because of the much smaller thermal velocity in 
the case of Calcium, the shear broadening is much stronger than for Hydrogen lines. This results in a stronger signature of spiral waves 
in the emission lines. 



phase. In the velocity coordinate space of Doppler tomo- 
grams, axisymmetric discs wiU produce axisymmetric cir- 
cular images, centered on the radial velocity of the white 
dwarf {Vx — 0,Vy = —Ki, with Ki the projected orbital ve- 
locity of the white dwarf). For each radius, a corresponding 
ring at the local disc velocity is produced, with the smallest 
velocities corresponding to the outer edge. 



3.2 Model I: High Mach number disc 

In Fig. ti we show the model predictions for the high Mach 
number CV disc of model I, viewed at an inclination of 80°. 
The spatial distribution of the line emission is plotted in 
the top left panel and displays a rather small disc, varying 
in radius between rout (0) — 0.55 — 0.65 tl^ , with two tightly 
wound spiral shock arms covering the disc down to small 
radii. If we plot the same distribution in velocity coordinates 
(top right), the typical ring shape appears, with the two 
shock arms visible as sharp arcs of enhanced line emissivity 
in the top and bottom areas. The outer disc edge provides a 
cut off of any emission with velocities lower than ~ 600 km 
s~^ while at high velocities, the inner disc, the tomogram 
becomes more symmetrical. The emission line profiles in the 
corresponding trailed spectrograms (Fig. S, bottom panels) 
feature two sharp peaks, reflecting the velocities of the outer 
disc gas. Two cases are plotted to illustrate the effect of shear 
broadening due to the finite thickness of the emission line 
region (Home & Marsh 1986). In the bottom left panels, 
the emission line layer is assumed to be very thin {AZ/H ~ 
0), so that shear broadening effects vanish and sharp and 



narrow double peaks are formed. This is very similar to the 
optical thin case. On the right is the case of maximum shear 
broadening, assuming the emission lines are formed across 
the full disc thickness such that AZ — H in eq. (5). Broad 
peaks with a V-shaped, rather than U-shaped valley between 
the peaks are formed. 

The presence of spiral arms slightly distorts the ve- 
locity of the peaks as a function of binary phase, vary- 
ing between Wpeak(0) = 690 - 750 km s~^ with respect to 
the white dwarf. On average, the double peak separation is 
(I'pcak) = 725 km s^^. This corresponds to a disc radius of 
Tout — G'Afi/(iipcak)^ = 0.61 r-Li , assumiug Keplerian veloci- 
ties. While the disc is slightly tidally distorted and elongated 
along the y-direction, the assumption of Keplerian veloci- 
ties still provides a reliable estimate of the average disc size 
in this case. 

Shear broadening can be particularly important for 
heavy elements. Their low thermal velocity results in a 
stronger local enhancement of the line flux due to velocity 
gradients. To illustrate this. Figure H shows the line profiles 
of a light element (A=l) versus that of Calcium (A=40) at 
an inclination of 45° . While for hydrogen like lines the shear 
broadening is no longer important due to foreshortening, the 
Calcium line displays strong anisotropic components origi- 
nating along the spiral arms. This increases the double peak 
variation from 8% in the case of no shear to 13% for maximal 
shear in heavy lines such as Calcium. 

The orbital variation in the separation of the double 
peaks, is the predominant signature that the disc is asym- 
metric and that the emissivity is not uniform. In the case 
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of considerable shear broadening, the line emission from the 
spiral shocks is additionally non-isotropic resulting in a bi- 
nary phase dependent variation of the peak strengths with a 
maximum at phases 0.25 and 0.75 as is visible in the emission 
line flux light curves to the right of the spectrograms. Such 
non-isotropic line emission cannot be modelled by Doppler 
tomography, since only the average line flux at a given ve- 
locity is provided, but can be identified in the emission line 
profiles directly as a signature of shocks in the outer disc. 
While saturated emission lines from heavy elements are thus 
good indicators of large velocity gradients, they are unfortu- 
nately usually weak lines and do not provide sufficient signal 
to noise for Doppler imaging in existing data sets. 
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3.3 Model II: dwarf nova outbursts 

By increasing the parameter of the small scale viscosity to 
a = 0.3 in Model II, viscous spreading becomes more effi- 
cient in pushing the accretion disc towards larger radii of 
rout(0) — 0.6 — 0.8 rLi (see Fig. 0, top left panel). The in- 
creased local viscous heating additionally produces higher 
disc temperatures, and thus lower Mach numbers, now of 
the order M ~ 5 — 20. This is reflected in the open geome- 
try of the two spiral shock arms in the outer disc. In Fig. 
we present the calculated emission line properties for model 
II, again at an inclination of i = 80°. 

It is difficult to assign a meaningful disc radius as the 
line emissivity is predominantly localized at the two spiral 
arms. The strong m=2 spiral wave pattern in the outer disc 
dominates the correspondiiig Doppler tomograms, shown in 
the top right panel of Fig. y. The line emission extends to- 
wards lower velocities {^^ 500 km s~^), partly because of the 
increased disc size, and partly because sub-Keplerian mo- 
tions are significant in these outer regions. Figure a shows 
the deviation from Keplerian velocities at a few representa- 
tive radii. While the disc is close to Keplerian in its inner 
regions (high velocities), the departures at the outer disc 
near v ~ 500 km s~^ are significant and of the order of 
-lOOkms"^ 

The line profiles are easily distinguishable from an ax- 
isymmetric disc. The m = 2 spiral pattern results in con- 
verging emission line peaks, with a cross over near phases 
0.25 and 0.75. This contrasts to the axisymmetric case, with 
a constant double peak separation, independent of orbital 
phase. Though a non-circular disc could also produce a vari- 
ation in the double peak separation with phase, it would 
not produce the existence of multiple peaks as one crosses 
from one spiral arm to the next, with a corresponding sud- 
den jump in velocity of several hundred km s^^, a unique 
signature of a spiral shock pattern. Again the large veloc- 
ity gradients are important at high inclinations where shear 
broadening can dominate. This produces a variation in the 
emission line light curve (Fig. bottom right panel) which 
is considerably larger compared to Model I. 



RECONSTRUCTING SPIRAL STRUCTURES 
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Figure 5. Non-Keplerian velocities in model II. The circles de- 
note the Kepler velocity at a few radii in the disc. The small ticks 
point to the actual disc velocity at that radius and azimuth. The 
tidal distortions become prominent in the outer disc, i.e. at low 
velocities in the inner regions of the Doppler maps. 



4.1 Realistic data quality 

In the previous section we discussed our model predictions 
limited by the grid resolution of our hydrodynamical models. 
We now include instrumental noise and finite spectral/time 
resolution in our predicted data to provide a realistic com- 
parison with observed data and constrain the necessary data 
requirements to reconstruct tidal arms in CV discs. 

Three instrumental effects play a role. First the wave- 
length resolution of the individual spectra, setting a limit to 
the velocities that can be resolved. The emission lines typ- 
ically span several thousands of km s~^, which need to be 
sampled by ample number of pixels. On the other hand, res- 
olutions much better than the width of the local line profile 
(~ 10 km s~^) are not providing any additional information. 
A resolution of ~ 20-80 km s~^ is therefor usually optimal, 
though is practically limited by signal to noise requirements. 

Secondly, the phase or time resolution has to be short 
enough to provide sufficient phase sampling across the bi- 
nary orbit to avoid artifacts in the image reconstruction (see 
also Marsh & Home 1988). The higher the resolution of the 
input data, the more projections are required to provide 
sufficient sampling and exploit the available resolution. To 
avoid under sampling the image one typically needs at least 
50 orbital phases. 

Finally, the signal to noise of the individual spectra will 
determine how well the data will constrain the image struc- 
ture. In the case of maximum entropy reconstruction as is 
the case here, the reconstructed image is a balance between 
fitting the data using maximum likelihood statistics and se- 
lecting the simplest image using the maximum entropy cri- 
terion. 

To test the importance of above mentioned effects, we 
start by calculating ideal line profiles (such as presented 
in Section pi) for an inclination of 45° at 50 orbital phases 
equally spaced over the binary orbit. The data is then con- 
volved with a Gaussian with FWHM equal to the desired 
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Figure 4. As Fig. ^ but now for tlie low Mach number disc of model II. 



Lightcurve 




instrumental resolution and binned in wavelength such that 
2 data pixels cover the instrumental profile. Poisson noise 
was subsequently introduced to achieve diflerent signal to 
noise levels in the line profiles as follows. We added a con- 
tinuum to our line profiles such that the line strength (line 
- continuum) equals about two times the continuum level, 
a typical value for the strong emission lines in dwarf novae. 
Poisson noise was then added to achieve the desired signal 



to noise level in the continuum. Using the synthetic noisy 
data, we construct maximum entropy Doppler tomograms 
using the same code that is used for the observations |j 



DOPPLER package developed bv Tom Ma rsh, see 
tittp : //www . astro . soton . ac . uk/~trm/sof tware . html 
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Figure 6. Tomograms and trailing spectra of the tiigii Macti number simulation (model I) for a strong emission line. From left to right 
the 2-pixel wavelength/velocity resolution is (40, 80, 80) km s~^ and the ratio of signal to noise is (50, 50, 15) in the continuum. 



4.2 High M model 

Figure ra shows the image reconstructions based on model 
I, at an inclination of 45° and including maximum shear 
broadening. The top panels show the simulated data, the 
left has a 2-pixel resolution of 40 km s~^, the middle and 
right have a resolution of 80 km s~^. The signal to noise in 
the left and middle panels is ~ 50 in the assumed contin- 
uum, while for the right panel it is degraded to ~ 15. Below 
each spectrogram the reconstructed Doppler tomograms are 
shown, all on the same scale. 

It is difficult to reconstruct a clear signature of the 
tightly wound spiral waves. Except for an enhancement of 
the disc emission in the top left and lower right, such discs 
will look very similar to axisymmetric discs, even if high 
quality data is available. As data quality is degraded the 
slight variation of the double peaks with binary phase be- 
comes difficult to measure. We conclude that we do not ex- 
pect to see obvious evidence of spiral waves in Doppler maps 
of small, low a accretion discs even if tightly wound spiral 
waves are present. 



4.3 Low M model 

Figure m plots similar reconstructions now based on model 
II. Resolution and signal to noise levels are identical to the 
previous model for comparison. The strong two armed spiral 
pattern is readily recovered even at low signal to noise and 
medium resolution (Fig. H, right). In these large hot discs, 
the tidal arms dominate the emission lines and are easily 
visible both in the line profiles directly and in Doppler to- 
mograms. 

To compare the two models, we show the asymmetric 
part of the reconstructed tomograms in Fig. 0. These are 
obtained by subtracting the median of the image at each ve- 
locity, measured from the white dwarf. The low Mach num- 
ber case (bottom) shows prominent spiral arms, whereas for 
model I, only a slight asymmetry is observed. 



4.4 Instrumental requirements 

We reconstructed disc images for various signal to noise lev- 
els, amounts of shear broadening, orbital inclinations and 
number of binary phases using Model II data. Based on these 
reconstructions, we found a signal to noise of ~ 15 in the line 
flux and a two pixel velocity resolution of ~80 km s~^ to be 
sufficient to convincingly reproduce the two armed pattern 
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on the disc. As we mentioned earlier, to avoid artifacts ~ 50 
orbital phases are desired. The instrumental requirements 
for reconstructing such open spirals are therefor fortunately 
relatively modest and suggest a search for such structure in 
CV discs is feasible and will provide strong constraints on 
the presence of spiral waves in these discs. 

Figure H illustrates some of the reconstructions ob- 
tained. The left two images illustrate the effect of shear 
broadening for high inclination cases. The disc image is sim- 
ilarly broadened in the case of strong shear broadening but 
still carries the two armed pattern. Middle panels illustrate 
the effect of a strong versus weak emission line. This basi- 
cally lowers the signal to noise of the data and affects the 
disc image in the same way. Finally the right panels show 
the strong artifacts due to poor phase sampling (top) or the 
lack of image structure due to poor resolution. 



4.5 Other emission line sources 

Our reconstructed images are those of the accretion disc 
contribution only. The irradiated companion star, the gas 
stream and the hot spot are well known other contributors 
to the emission line flux. In the case of the companion star, 
its contribution is confined to a Roche lobe shape centered at 
the radial velocity of the companion, at velocities where no 



disc emission is expected. More demanding are the stream- 
disc interactions, which can contribute to a wider range of 
velocities overlapping the accretion disc, though the com- 
pact impact region is usually well defined. Complicated are 
disc outfiows, which can be responsible for the single peaked 
emission lines observed in some systems. These issues have 
to be born in mind when interpreting emission line data sets 
(e.g. Steeghs, Home, Marsh & Donati 1996). 



5 DISCUSSION 

The Doppler images and line profiles presented in this paper, 
confirm that the observed disc structure in tomograms of IP 
Pegasi can be reproduced by tidally excited spiral arms in 
the disc. The presence of a dominating two armed structure 
in early outburst, when the disc is hot and large, compared 
to a fairly symmetric disc in quiescence matches extremely 
well 'with the model predictions presented here. Armitage 

6 Murray (1998) reach a similar conclusion based on 3-D 
SPH simulations of the disc fiow in IP Pegasi. A spiral den- 
sity pattern, matching the kinematics of the observations, is 
formed when a is increased and the disc expands. However, 
these arms appear to be transient in their simulations on a 
viscous time scale. In the models of Stehle (1999) the spiral 
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Figure 9. Some example reconstructions for Model II data. Left is the case of high inclination (80°) without (top) and with maximal 
shear broadening. Middle a reconstruction of a strong emission line (twice the continuum) with a continuum signal to noise of 30 (top) 
compared to a weak line (15% above continuum) at the same S/N. Top right shows artifacts introduced by poor phase sampling, in this 
case 10 bins even at high signal to noise (50). Bottom right panel is finally the case of insufficient resolution (150 km s~^) to reconstruct 
the disc structure. 



pattern is stable in the co-rotating frame, though its prop- 
erties will change as the outburst progresses and especially 
when the cooling wave starts to turn the disc into quies- 
cence again. Such time dependent behaviour can be tested 
by Doppler mapping of the disc as the outburst progresses. 
Godon, Livio & Lubow (1998) presented 2-D disc simula- 
tions using Fourier-Chebyshev spectral methods. They at- 
tempt to reproduce the observations of IP Pegasi. They find, 
however, that the spiral pattern resembles the observations 
only for very high temperatures. Unfortunately their a val- 
ues were rather low {a = 0.1 at most) and thus viscous 
spreading was less efficient compared to our simulations. As 
a consequence their discs where smaller by up to ~ 50% 
compared to our model II simulation. Indeed, our model I 
simulation shows that in the case of low a values little evi- 
dence of spiral arms in the emission lines is expected. How- 
ever, with efficient viscous spreading the disc becomes large 
enough to result in a spiral pattern comparable to observed 
Doppler tomograms. 

The behaviour of the emission lines in the presence 
of strong spiral waves is very different from axisymmetric 
or elliptically distorted accretion discs without spiral shock 
waves. In axisymmetric discs we expect a constant separa- 
tion between the peaks, and in elliptical discs a smooth, sinu- 



soidal variation with orbital phase, however without sudden 
jumps in the peak separation or the presence of multiple 
peaks arising form successive spiral arms such as in Fig. 
m Even though spiral arms in close binary accretion discs 
have been predicted more than 10 years ago (Sawada et al. 
1986), and Doppler tomography provides a way of resolving 
the dynamics of accretion discs, a clear detection of a spiral 
pattern was only achieved in a recent outburst of IP Pegasi. 
Although confirmed since in this particular object (Harlaftis 
et al. 1999), detection of such a pattern in other systems is 
crucial. 

We showed, however, that in small, cool discs, the ob- 
servational signatures left by tightly wound spiral waves are 
small. Indeed, as Figure H shows, such discs, will closely re- 
semble axisymmetric discs. The data on IP Pegasi and our 
simulations suggest that the increase in disc size and tem- 
perature during dwarf novae outbursts results in a promi- 
nent spiral pattern in the disc that can be recovered using 
tomography. The number of Doppler images of dwarf novae 
in outburst is unfortunately small. Marsh & Home (1990) 
provide Doppler images of IP Pegasi in both quiescence and 
outburst. In particular their outburst map of Hell (A4686) 
four days into outburst show a two armed enhancement in 
the disc at the correct positions in the tomogram. Their res- 
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Figure 8. The asymmetric parts of the reconstructed imaged for 
the two models. In the high M case, only a slight asymmetry is 
produced in the top left and lower right regions. In the low M 
case, the image is dominated by the two armed spiral structure. 



olution was unfortunately rather low (~130 km s~^). Piche 
& Szkody (1989) note that their outburst data of IP Peg is 
not compatible with a symmetric Keplerian disc, but no disc 
image is available. Similarly, Hessman (1989) reports mea- 
surements of the emission line peak velocities of IP Pegasi 
at the end of an outburst which shows a phase dependence 
that suggests a strong asymmetry in the outer disc persists 
well into the outburst, though is much weaker than in the 
initial phases. Outburst tomograms of SS Cyg (Steeghs et al. 
1996) show a highly non-axisymmetric disc, again with en- 
hanced emission in the top right and lower left areas. On the 
other hand, tomography of the disc in OY Car in outburst 
(Harlaftis & Marsh 1996) shows some evidence for the gas 
stream and its impact on the disc, but no clear spiral disc 
structure. Whether this is related to the fact that OY Car is 
a SU UMa type dwarf nova, will be explored in the future. 



Nova-like variables provide permanent high mass transfer 
accretion discs, and we would therefore expect tidal struc- 
tures to be important. Unfortunately the emission lines from 
these systems are poorly understood. The lines are single 
peaked and contain S-wave components of unknown origin. 
The exception could be V347 Pup, the only eclipsing nova- 
like variable with clear double peaked emission line profiles 
originating in the disc (Still et al. 1998, Steeghs et al., in 
preparation). The disc appears to be tidally distorted, and 
more data may confirm the presence of spiral structure in 
this system. 

Conclusive remarks on the presence of spiral waves in 
accretion discs are therefor difficult to make with the sparse 
amount of outburst data available at this time. We have 
discussed the properties of CV discs since these are still the 
best laboratories to study the properties of thin, a type 
accretion discs where most heating is due to viscous dissi- 
pation. In particular, high quality spectroscopic studies of 
dwarf novae throughout their outburst cycles are invaluable 
in understanding the physics of thin accretion discs. It is 
important to confirm with further observations that the spi- 
ral density distribution in close binary accretion discs is a 
common phenomenon and during which source states they 
can be observed. 

The significance of tidal waves in the discs of X-ray bi- 
naries (Blondin & Owen 1997, Owen & Blondin 1997, Lan- 
zafame & Belvedere 1998) is also considerable, since the 
Mach numbers in these discs are much lower and the rel- 
evance of density waves to the total angular momentum 
budget of the disc is therefore greater compared to CVs. In 
due course this should be accessible to similar observational 
tests. 



6 SUMMARY 

We studied the emission line profiles of close binary accretion 
discs in cases where the disc pattern is dominated by promi- 
nent spiral shock arms excited by the tidal forces of the sec- 
ondary star. As underlying disc models we used the grid of 
full hydrodynamic CV disc calculations of Stehle (1999). We 
model the emission lines using the theory of Home (1995), 
including the effect of shear broadening for saturated lines. 
Doppler tomograms and trailing spectra where constructed 
and instrumental effects were included to assess the data 
quality required to reconstruct spiral structure in CV discs. 

We found that for high Mach number discs (M = 
v^/ca ~15 - 30) the spiral shock arms are too tightly wound 
to leave significant fingerprints in the emission line structure 
of these small discs. As a result the peak separation of the 
double peaked emission line profile varies only by ~ 8% with 
binary phase. We therefore do not expect to detect strong 
spiral arms in quiescent CV discs. 

For low Mach number discs (M ~ 5 - 20), i.e. for CV 
accretion discs in outburst, the line emissivity is dominated 
by the spiral shock arms. The two armed spiral pattern is 
identified by the presence of converging peaks in the line 
profiles. This signature in the trailing spectrum is expected 
to be sufficiently prominent that a signal to noise of ~ 15 
per pixel, a wavelength resolution of ~ 80 km s~^ and a 
time resolution of 50 spectra per binary orbit is sufficient to 
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verify the presence of spiral shock arms in binary accretion 
discs through Doppler tomography of strong emission lines. 

Our simulations successfully reproduce the emission 
line properties of the disc structure observed in IP Pegasi. 
The detection of strong two armed spiral arms in outburst 
(Steeghs et al. 1997), but no clear disc structure in quies- 
cence (Marsh & Home 1990), is convincingly reproduced by 
tidal spiral waves in the accretion disc. 

The prospect of a new way to probe the local physics of 
discs through the dynamics of these tidal spirals demands 
more observational efforts on the detailed properties of high- 
mass transfer accretion discs through indirect imaging. 
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